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We report on the mixed-orbital ground state compounds analogous to the mixed-valence ones. 
Orbital doublet and singlet ground states of the Fe 2+ ion displayed in Mossbauer spectra of the 
Prussian blue analogues Aj / Fe3_ H [Co(CN)e]2 -ccEbO (y ^ 1) are interconverting to each other as 
temperature changes for A = K, Na. Relative weight of orbital singlet ground states increases with 
lowering temperature. In the alcali-free cobalticyanides, the Mossbauer spectra are dominated by 
quadrupole-splt doublets of the same origin with a large splitting, characteristic of the Fe 2+ species 
coordinated by two (at least) oxygen atoms of water molecules in cis-configuration. Single-type 
octahedral Fe 2+ coordination is isolated in A-filled cobalticyanides for A=Rb, Cs to be characterized 
by the narrower quadrupole spectrum associated with the orbital doublet ground states. Ionic 
exchange of K for Cs in KFe[Co(CN)e] -H^O results in predominance of the orbital singlet, i.e., in 
enhanced anisotropy of the Fe 2+ valence electrons. We found hence that the Fe 2+ charge distribution 
can be modified together with crystal size at the step of synthesis owing to distortion isomerism. 

PACS numbers: 71.70.-d Level splitting and interactions; 76.80.+y Mossbauer effect; 



I. INTRODUCTION 

Materials enclosing the transition metal elements are 
characterized by a wide diversity of ground states 
whether this be the quantum states of collective electron 
systems or local states. In latter case, the transitions be- 
tween disparate ground states of 3d ions can be driven by 
changes in temperature, chemical doping, or by external 
stimuli. Phase transitions are thus induced that involve 
the spin, charge, lattice and orbital degrees of freedom. 
Charge transfer, spin crossover and linkage isomerism are 
the three types of transitions requiring a complex elucida- 
tion of the roles of correspondingdegrees of freedom to be 
distinguished from each othepl"2. The long-range struc- 
tural aspects display their role via succession of transi- 
tion steps as they evidence the interconversions between 
various local states mediated by elastic interactions^^. 
Charge transfer from one site to another could be a pre- 
requisite of the concomitant spin-state transitiorPpl Or- 
bital degree of freedom is responsible for a long-range- 
distortion symmetry lowering, provided that orbital or- 
dering is coupled with the mixed-valence charge ordering 
transitiorP. Orbital ground state reversaPp^ is an exam- 
ple of yet another type of transitions strongly coupled to 
the lattice and, possibly, to some other degrees of free- 
dom. Simultaneous presence of the orbital singlet and 
orbital doublet ground states, familiar in case of Fe 2+ 
iorP3] could be an indication of proximity of such a tran- 
sition. 

Prussian blue (PB) analogues constitute a family of 
simple cubic compounds which excited a great interest 
recently owing to the family members exemplifying the 
spin crossover, charge transfer and linkage isomerisrrPP 
Here, in the same family, Aj,Fe3_j / [Co(CN)e] 2 -nH^O 
(y ^ 1, A = K, Na, Rb, Cs), we report on mixing and 
intercoverting the orbital ground states of Fe 2+ . 

Metal-organic framework of the PB analogues is 



formed by the 3d metal cations and the cyanometal- 
late complex anions (of the [Co(CN)g] 3_ -type) arranged 
into a NaCl-like 3D checkerboard. In major group of 
these compounds, the valence and spin differ between 
the states of metal ions occupying the cationic node and 
the anion core. The framework complex anions are typi- 
cally having higher valences than the cations, in this case, 
the framework is bearing a net negative charge, which is 
balanced either by insertion of alcali (A) cations into the 
framework interstitials or by spontaneous generation of 
the large-size vacancies on the anionic nodes. Crystalliza- 
tion and zeolitic water molecules (C-water and Z-water) 
fill the cavities at the vacant sites of complex anions and 
A-cations, respectively. 

The local electronic configurations are probed by the 
quadrupole interactions of a 57 Fe Mossbauer nuclei sens- 
ing the electric field gradient from valent electrons and 
lattice charges surrounding the Fe 2+ ion. Mossbauer 
spectra of ferrous hexacyanocobaltates were reported 
previously only for the alcali-free phases (y = 0]f^p^. 
The spectra can be fitted by two or more quadrupole 
doublets implying the existence of a number of differ- 
ent states of the high-spin Fe 2+ ions. Paired doublets, 
both originating from the cationic sites outside the low- 
spin complex [M(CN)e] were also observed in Mossbauer 
spectra of some other ferrous hexacyanometallates, e.g., 
for M=Cr (y = 0)P. The early Mossbauer studjD^ has 
proposed that the water molecules in cis- and trans- con- 
figurations could contribute to the appearance of multiple 
sites for Fe 2+ . In case of M=Cr, the spin-state transition 
of Fe 2+ was observed only for y ^ (A=Cs)P^. It was 
conjectured^ that filling the structure with A ions in- 
duces a stronger ligand field on the Fe 2+ ion that causes 
the spin transition. In this study of the local electronic 
configurations in AFe[Co(CN) 6 ]-H 2 ( A = K, Na, Rb, 
Cs) and related alcali-free form we show that the ligand 
field stabilizes the singular ground state (orbital doublet) 
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only for A=Cs and Rb. In the cases of A=Na, K, and 
y = we report on the coexisting orbital singlet and 
orbital doublet ground states for Fe 2+ and on the inter- 
conversions between them induced by changes of temper- 
ature. 



[a) F-43m 
Alcali-filled ordered 



(b) Pm-3m 
Alcalr-free ordered 



II. 



EXPERIMENTAL 



RbFe[Co(CN) 6 ]-H 2 and CsFe[Co(CN) 6 ]-H 2 were 
prepared mixing the 0.025 M aqueous solution of 
K 3 Co(CN) 6 with 0.025 M solutions of FeCl 2 -4H 2 in 
excess of RbCl or CsCl. Similarly synthesized sam- 
ples in excess of KI and NaCl gave the compounds 
K . 8 Fei.i[Co(CN) 6 ]-H 2 Oand Na .7Fe 1 . 1 5[Co(CN) 6 ]-H 2 O, 
respectively. Analytic methods of inductively-coupled 
plasma and x-ray-fluorescent analysis gave the results 
consistent with the site occupancies obtained from Ri- 
etveld analysis. The same compounds with different 
crystal size were synthesized for A=Rb and Cs via 
the ionic exchange using K .8Fei.i[Co(CN) 6 ]-H 2 O as a 
sorbent (Table I). Transformation of this sorbent into 
CsFe[Co(CN) 6 ] was detected by Rietveld analysis and 
confirmed by analytic techniques. 

Table I. Lattice parameters in the precipitated and 
ion-exchanged hexacyanocobaltates at room tempera- 
ture. 



Compound 


Parameter*, nm 


Na . 7 Fe 1 . 15 [Co(CN) 6 ]-H 2 O 

K . 8 Fe 1 . 1 [Co(CN) 6 ]-H 2 O 

KFe[Co(CN) 6 ]-H 2 

RbFe[Co(CN) 6 ]-H 2 

CsFe[Co(CN) 6 ]-H 2 

Fe[Co(CN) 6 ] . 67 • 5H 2 

Cs-exchanged (b) at 20°C 


1.03278(8) 
1.03115(3) 
1.03107(3) 
1.03455(7) 
1.03682(8) 
1.02778(8) 
1.03628(8) 



*Refmed according to Fig.l (c) and (d) structure mod- 
els using the Fm3m space group. 

X-ray Diffraction (XRD) patterns were measured on 
a PW3040/60 X'pert PRO (PANalytical) diffractometer 
using Cu Ka source. The atomic ratio Fe/Co and K/Co 
were determined by inductively coupled plasma atomic 
emission spectrometer (ICP-AES) on an IRIS Interpid II 
XSP Instrument (Thermo Electron Corporation). The 
x-ray fluorescent analysis was carried out using the XRF 
Axios X (PANalytical) spectrometer. The water content 
in the samples was estimated by thermogravimetric anal- 
ysis carried out using a Setaram Setsys 16/18 analyser 
measuring the weight loss in an air flow of 25 ml/min 
with a heating rate of 10 K/min. Magnetization of the 
paramagnetic compounds was measured with a SQUID 
magnetometer in the range between 2 K to 350K. 

Mossbauer spectra were measured with a Topologic 
500 A spectrometer at the temperatures between 78 K 
and room temperature. The measurements at heating 
the sample above room temperature were performed in 
a Wissel GMBH furnace. Isomer shifts (IS) are reported 
relative to a-Fe at room temperature. 




(c) Fm-3m 
Alcali-filled disordered 



(d) Fm-3m 
Alcali-free disordered 
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FIG. 1: Crystal structures and symmetry groups of the four 
modifications of Prussian blue analogues: alcali-filled ordered 
F43m (a), alcali-free ordered Pm3m (b), alcali filled disor- 
dered Fm3m (c) and alcali-free disordered Fm3m (d). 



III. RESULTS AND DISCUSSION 
A. Structural considerations 

The A-filled and A-free Prussian Blues are two types 
of the cubic frameworks, described by structure models 
of Keggin and Miles^ and Ludf^, respectively. Both of 
them exhibit the variable degrees of randomness in the 
distribution of the A ions over the framework interstitials 
and in the distribution of the water-filled anionic vacan- 
cies. In both cases, the random structures are described 
by the symmetry group Fm3m. When the A-cations and 
Z-waters occupy the interstitials in an ordered fashion 
the symmetry turns into F43m. The F43m symmetry ap- 
pears, for example, in CsFe[Cr(CN)6]-H 2 15 . It means 
in the ideal case that the system of framework intersti- 
tials splits into two sublattices (1:1), as shown in Fig 
1(a). In reality, each of these (1:1) sublattices is not rep- 
resented by solely A-ions and Z-waters, but occupied just 
unequally with these two species. 

In the A-free structure, another sublattice may split 
into two subsystems that is the anionic sublattice. When 
one proceeds from the disordered to ordered placement 



of anionic vacancies the symmetry turns from Fm3m to 
Pm3m giving birth to two dissimilar sublattices 18 . In the 
Pm3m structure, two born sites are abundant with the 
ratio of 3:1. Four relevant lattice cells are shown in Fig.l 
for the pairs of A-free and A-filled structures. 

Rietvcld refinements were carried out using the above 
mentioned four structural models and three symmetry 
groups known from the literature^HII]. These treatments 
have shown that the most symmetric Fm3m model is 
sufficient to fit all the observed reflections. When fitted 
with lower symmetries the XRD least-square analysis has 
resulted in stagnant goodness of fit, unimproved, despite 
larger number of fitting parameters. The intensity ratio 
^200/^220 changes dramatically when K, Rb or Cs ions 
enter into the structure, as reported previously®. The 
change becomes very large as the number of electrons in 
an alcali ion increases, allowing for Cs a rough estimation 
of its content even prior to Rietveld analysis. 

The refined lattice parameters at ambient tempera- 
ture are listed in Table I. In all compounds, the lat- 
tice parameter is close to that of the low-spin phase of 
CsFe[Cr(CN) 6 ]-H 2 CPl, however, it is found fr om mag- 
netic measurements that the addressed in Table I com- 
pounds enclose the Fe 2+ ions in the high-spin state that 
is unchanged down to 2 K. 



B. Doublet and singlet orbital ground states of 
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Fe in PB analogues 



FIG. 2: Formation of the orbital doublet ground state of Fe 2+ 
in the "5+1" environment. In a Prussian Blue analogue, the 
valence spin-down electron clouds are depicted for the Fe 2+ 
species coordinated by 5 nitrogens and one water molecule. 
It is shown that two orbitals of the t2 g group (d xz ,d yz ) retain 
the equivalent energy (a,b), whereas the third one (c) splits 
off upwards. A single spin-down electron occupies the double- 
degenerate level forming the doublet orbital ground state. 



1. Splitting the d-orbital energy level in anisotropic 
environment 

In the fully A-filled PB analogues, the Fe 2+ ions are co- 
ordinated with six nitrogens, however, at partial filling at 
least one C-water molecule enters the first coordination 
sphere of Fe 2+ . From the combined scheme of compen- 
sating the charge misbalance between Fe 2+ cations and 
Co(CN)g] 3 " anions this is evidently the case of our sam- 
ples (Table I). Electronic configuration of a high-spin 
d 6 ion implies only one electron extraneous with respect 
to isotropic half-closed d-shell composed of five spins-up 
electrons. In an octahedral environment, the group of 
three t2 g orbitals is available for this sixth electron. We 
show in Fig. 2 the orientations of these orbitals relative 
to the water molecule. 

The values of chemical shifts for both doublets in Fig. 
3 are in the range 1.1 mm/s to 1.15 mm/s allowing us 
to identify indeed both doublets with the high-spin Fe 2+ . 
In case of high-spin Fe 2+ , the valence contribution of the 
electric field gradient (EFG) is predominant and exceeds 
the lattice EFG term by an order of mag nitudeP. The 
t2 g group is thus the source of valence term of EFG that 
is the main contribution of the observed quadrupole split- 
ting. 

One of the orbitals of the t2 g group in Fig. 2 (a,b,c) 
is perpendicular to the direction pointing towards the 
water molecule. Corresponding energy level would not 



be affected by the ligand replacement. The lobes of two 
other orbitals (d xz and d yz in Fig. 2) are at 45° with re- 
spect to oxygen direction. Owing to reduced repulsion 
between oxygen and d xz or d yz electrons corresponding 
levels would have lower energy. In the ground state of 
Fe 2+ ion, these d xz and d yz share the sixth d-electron 
(spin-down electron), while five spin- up d-electrons span 
isotropically. The third {d xy ) orbital would split off its 
energy level upwards from the ground state, formed by 
two degenerate d xz and d yz orbitals as shown in Fig. 
2(d). 

2. Temperature dependence of quadrupole splitting in 
A-filled cobalticyanides (A=Na,K): 

We observe in Fig. 3 two doublets, both of them 
having the temperature-dependent quadrupole splitting. 
In addition, the ratio of doublet areas is also strongly 
temperature-dependent. The external doublet is domi- 
nating at low temperatures. As the temperature is raised 
an inverted doublet area ratio is observed. Mossbauer 
spectra show that the triple degeneracy of t 2g is lifted 
over either in favor of the orbital singlet or in favor 
of orbital doublet ground states. The charge distribu- 
tion is less anisotropic when the spin-down electron is 
shared between two orbitals as shown in Fig. 2. We as- 
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(a)Na 07 Fe 115 [Co(CN) 6 ]xH 2 O (b) K^Fe, ^[CofCNy-xHjO 




o Temperature (K) o Temperature (K) 



FIG. 3: Mossbauer spectra of the A- filled ferrous hexa- 
cyanocobaltates for A=Na and K in temperature range be- 
tween 80 K and 280 K. In the bottom panels, the experi- 
mental temperature dependences of the quadrupole splittings 
determined from least-squares fitting the Mossbauer spectra 
are shown by triangle symbols. The solid lines represents 
AEq(T) fitted using Eq.(l) with the parameters Ai and A2 
indicated. 



cribe such a charge state to the internal doublet of Fig. 3. 
Larger charge anisotropy originates from localizing the 
sixth electron of Fe 2+ on a single d-orbital (external QS 
doublet). From the QS ratio the valence EFG term is 
expected to be twice larger for the orbital singlet than 
for the orbital doublet. The literature overview con- 
firms indeed that in Fe 2+ the orbital singlet produces 
twice larger qua drupol e splitting than the orbital dou- 
blet ground statfPSEDEI 

In the bottom of Fig. 3, we fitted the values of 
AEq derived from Mossbauer spectra with the theoreti- 
cal expression^ 



AE Q (T) = 



-2Ax -2A 2 ~Aj -A 2 -(A!+A 2 ) 

AE Q EaT =S£ 

1 + e k B T + e k B T 

(1) 

Here AE'q is the quadrupole splitting at zero tempera- 
ture, which can be as high as 3.7 mm/s for Fe 2+ ^. The 
fitted values of the energy gaps Ai, A2 are also shown in 
Fig.3. 

The ratio of the fitted values Ai : A2 is small (~ 0) and 
large (~ I) for the internal and external QS subspectra, 
respectively. This is consistent with the subspectra as- 
signment to the orbital doublet and orbital singlet ground 
states identified according to their QS values. Very small 
splitting Ai ( <C A2 ), similar to fitting error, is allowed 
between the d xy and d xz levels. Even if this degener- 
acy is locally lifted by some remote disorder, such small 
splitting Ai remains unobservable until very low temper- 
atures. On the other hand, the effect of larger A2 is the 
robust steepness of the QS vs. T dependence troughout a 
broad temperature range. Similar splitting of 340 cm -1 
(ra 490 K) was estimated in Fe 3 [Co(CN) 6 ]2 ■ 5H 2 by 
Rasmussen and Meyers^. 

3. Classification scheme of the Fe 2+ sites by the symmetry 
of environment 

Our Rietveld refinements in both A-filled and A-free 
series of cobalticyanides have revealed that their struc- 
tures are best described by the symmetry Fm3m, allow- 
ing for the Fe 2+ ion a number of different coordination en- 
vironments. The Fe 2+ sites can be classified by the num- 
ber and configuration of N-ligands replaced with H2O. In 
these Fe 2+ sites, the lowest in energy state depends on 
the H2O inclusion in first coordination sphere, however, 
when H2O is not included, the ground state is deter- 
mined by the local distortion. Owing to so-called distor- 
tion isomerism^ the ground state becomes not uniquely 
defined. Shown in Fig. 4 are the cases of distortion iso- 
merism (a), orbital singlet (b), and orbital doublet (c). 

The distortion isomerism is the phenomenon closely 
related to the inability for the sixth (spin-down) electron 
to be distributed isotropically over the triply degener- 
ate t2 9 orbital. Such an overmuch triple degeneracy is 
universally lifted over by a local distortion of one or an- 
other type. In a FeN 6 octahedron, the compression along 
tetragonal axis would stabilize the orbital singlet lowest 
(a-I), while elongation along tetragonal axis stabilizes 
the orbital doublet as a ground state (a-2). Opposite 
effects are expected for the compression and elongation 
along the trigonal axis of the octahedron^!. Self-coupling 



5 




xy,xz 



FIG. 4: Energy level diagrams and corresponding coordina- 
tions for the orbital singlet (a-1) and (b) and orbital doublet 
(a-2) and (c) ground states of the Fe 2+ ion in ferrous hex- 
acyanocobaltate. In a-1 (a-2) the diagram of level splitting 
for the distorted FeNg ocahedron is shown for the elongation 
(compression) along tetragonal (trigonal) axis. 



between electron localization and local distortion elimi- 
nates the case of zero EFG valence term that is never 
observed. 

This physical picture is a key to understanding 
the temperature variation of the abundance of two 
Mossbauer doublets in Fig. 3. Assuming that the 6N- 
coordinated Fe 2+ sites are contributing to both doublets, 
we observe that the relative abundance of the orbital 
doublet and orbital singlet states is varying with tem- 
perature owing to varied lattice distortion. Matsuda et 
aP^ have reported that the PB lattice possess an intrin- 
sic instability with respect to rotations of the rigid units 
of the [Co(CN)6] 3- -type. The rotational displacements 
generate spontaneously to stabilize the excited charge- 
transferred states, which arises owing to strong lattice- 
electron coupling. The same rotational displacements 
cause also the negative thermal expansion^. Recently, 
the thermal expansion in Fe3[Co(CN) 6 ] 2 • nH 2 was re- 
ported to be indeed negative, at least above 200 Kp3. 

We observe thus that the orbital ground state reversal 
accompanies the distortion isomerism. One or another 
of coexisting distortion isomers prevails in a particular 
temperature range. Largest charge anisotropy occurs at 
lower temperatures owing to localizing the sixth electron 



of Fe 2+ on a single d-orbital. The charge distribution is 
less anisotropic at higher temperatures when the spin- 
down electron is shared between two orbitals. The tem- 
perature variation of the abundance of two Mossbauer 
doublets is very smooth. Similar temperature-varied lo- 
cal distortion concurrence that engenders the distortion 
isomerism was observed previously^. The glasslike tran- 
sition driven by supercooling AT was characterized^ by 
an equilibration time which decreases smoothly with tem- 
perature. The distortion switch-over is local and dis- 
tributed over a broad temperature range. 

With decreasing the content of A-ions the anionic va- 
cancies (water-filled pores) are generated for charge com- 
pensation. When the pore concentration is small only 
populated are the Fe 2+ sites which contain one H2O lig- 
and as shown in Figs. 2 and 4(c). Two pores in a vicinity 
of a Fe 2+ species occur either in trans- or in cis- config- 
urations. Consider in Fig. 4(c), for example, the ground 
state of the trans- Fe 2+ . If we denote by x the direction 
from the central Fe 2+ ion towards two trans- H2O ligands 
then the lobes of both xy and xz orbitals would be point- 
ing to directions having an angle of 45° with the x axis. 
The lobes of the third t^g orbital (yz) would be point- 
ing to the direction orthogonal to the x-axis. Therefore 
the ground state of the Fe 2+ ion in such FeN4(H20)2 
octahedron would be the orbital doublet. This is the 
same ground state as shown in Fig. 2 for the Fe 2+ ion 
in the FeN 5 H 2 octahedron with a single water in first 
coordination sphere. Orbital singlet in Fig. 4(b) is the 
ground state corresponding to the cis-configuration of 
FeN4(H20) 2 and to the octahedron with three in-plane 
waters. 

The H2O placement into the first coordination sphere 
of Fe 2+ suppresses the distortion isomerism and stabi- 
lizes a singular ground state of Fe 2+ . In our samples 
Nao. 7 Fe 1 . 1 5[Co(CN) 6 ]-H 2 and K . 8 Fe 1 . 1 [Co(CN) 6 ]-H 2 O 
the density of anionic vacancies is rather small and the 
population of the Fe 2+ ions surrounded by two vacancies 
is also small. In these samples the inner doublet comes 
from the Fe 2+ orbital doublet ground state associated 
with the FeNsH^O octahedron and with one of distor- 
tion isomers of FeNg. The outer doublet comes mainly 
from another distortion isomer of FeNg having the orbital 
singlet ground state. 



4- A-filled cobalticyamdes (A=Rb,Cs) 

Another manifestation of distortion isomerism is the 
difference between the Mossbauer spectra of the A- 
filled cobalticyanides obtained for Rb and Cs by two 
different procedures. Using the standard technique, 
the A-filled cobalticyanides were obtained by precipita- 
tion method. These compounds have demonstrated the 
single-component Mossbauer spectra (a) and (b) in Fig. 5, 
identified with orbital doublet ground state. On the other 
hand, when the method of ionic exchange was applied to 
replace K by Cs in KFe[Co(CN) 6 ]-H 2 0, the second QS 
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FIG. 5: Mossbauer spectra of the A-filled ferrous hexa- 
cyanocobaltates obtained in three different procedures: pre- 
cipitation (a-c,e), and ionic-exchange (d). The spectra fitting 
parameters are listed in Table II. 



subspectrum, identified with the orbital singlet, did not 
disappear. On the contrary, we observed the strong en- 
hancement of the outer doublet with A£ ! q=1.8 mm/s, 
from 32% to 64%. In addition, there appears a compo- 
nent (13%) with extremely large A£q of 3 mm/s. It 
is assigned to another orbital singlet state with larger 
energy level splitting parameters A 1; A 2 . The difference 
between two external doublets of the sample (d) is a man- 
ifestation another type of distortion isomerism unrelated 
to switching between different orbital ground states^. 

Rietveld refinements resulted in similar compositions 
for two different syntheses. The intensity of x-ray diffrac- 
tion has very similar patterns, except linewidth differ- 
ence. Larger-size crystallites of CsFe[Co(CN) 6 ]-H 2 are 
obtained via ionic exchange. In this method, the crystal 
size of is CsFe[Co(CN)e]-H 2 is determined by the size 
of the KFe[Co(CN)g]-H 2 precursor crystals. The target 
compound CsFe[Co(CN)e]-H20 with such large crystals 
could not be obtained by the standard method of precip- 
itation from a solute. Despite of similarity of the x-ray 
patterns, the QS of the ion-exchanged and precipitated 
samples of CsFe[Co(CN)e]-H 2 indicates that the valent 
spin-down electron has a very different distributions over 
the orbitals of t 2g group (Table II). 



Table II. Parameters of Mossbauer spectra (Fig. 5) 
in as-synthesized (a-c) and (e) and Cs-sorbed(d) hexa- 
cyanocobaltates: isomer shift (IS) S, quadrupole splitting 
(QS) AEq, linewidths r(mm/s), and spectral abundance 
(%). 



Compound 


5(mm/s) 




r(mm/s) 


% 






(mm/s) 






(a):RbFe[Co(CN) 6 ] 


1.091(5) 


0.992(8) 


0.60(1) 


100 


•H 2 










(b):CsFe[Co(CN) 6 ] 


1.116(3) 


1.077(4) 


0.71(1) 


100 


•H 2 










(c):KFe[Co(CN) 6 ] 


1.144(3) 


1.77(1) 


0.44(2) 


32 


TT f \ 

•H 2 


1.104(2) 


1.03(1) 


0.52(2) 


64 




0.39(1) 


0.66(2) 


0.23(3) 
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(d): ion-exchanged 


1.111(5) 


1.80(3) 


0.69(4) 


64 


K for Cs in (c) 


1.084(9) 


0.94(3) 


0.45(4) 


23 




1.123(8) 


2.98(2) 


0.33(4) 


13 


(e):Fe[Co(CN) 6 ] 2/3 


1.114(2) 


1.65(1) 


0.45(1) 


64 


•5H 2 


1.104(2) 


0.86(1) 


0.38(1) 


31 




0.33(1) 


0.23(2) 


0.26(3) 
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We observe in the standard (not ion-exchanged) com- 
pounds that only one type of the octahedron distortion 
is preferred when the large-size A-ions occupy the in- 
terstitials, and this type of distortion stabilizes the or- 
bital doublet lowest. On the other hand, in the case 
of A=K or A=Na the appearance of the outer doublet 
comes into play. In this series of PB analogues, sev- 
eral related aspects of structure and properties turn out 
to be A-size-dependent. Whereas in the case Cs(Rb)- 
containing PB analogue compounds the Cs(Rb) ions are 
situated in the centers of the interstitial cubes, in the case 
of some K(Na)-containing compounds the K(Na) ions are 
reported to be eccentrically located^. The mobility of al- 
kali metal ions in the channels of the hexacyanometalate 
structure decrease dramatically with increasing the alcali 
ionic radius^. The irreversible fixation of the less mobile 
ions ensures the Cs-selectivity for the cobalticyanide sor- 
bents, allowing the ionic exchange to be used in practice 
for 137 Cs removal from radioactive wastes 30 "^. 



5. Alcah-free ferrous cobalticyanide 

The isomer shifts in Mossbauer spectrum the A-free 
Fe[Co(CN) 6 ] 2 / 3 (Fig.5e) indicates predominantly (95%) 
the high-spin Fe 2+ . The ground states of the major part 
of the Fe 2+ species (64%) are the orbital singlet states. 
Third minor doublet is assigned to Fe 3+ according to its 
isomer shift and magnetic measurements. Refinements 
of the structure of Fei.5[Co(CN)6]-H 2 were useful to 
confirm the assignments of three Fe 2+ Mossbauer sub- 
spectra to the defined structural sites. Let us consider 
first the structure m odels of archetype Prussian Blue in- 
vestigated previous!} * 18 * 32 !. The concentration of pores in 
the A-free hexacyanocobaltate is very large, therefore, 
they cannot be considered as isolated ones. All kinds of 
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associations between the water-filled pores underlie the 
resulting abundance of the Fe 2+ species in various coor- 
dinations. Interactions between the pores at the stage of 
synthesis may result in various distributions of the pores. 
For the Prussian Blue Fe 4 [Fe(CN) 6 ] 3 • 14H 2 0, depending 
on the growth conditions, Buser et at^ could synthesize 
the single crystals with varied degree of disorder. Struc- 
tures of several crystals were refined using the Fm3m 
model, others were best refined using the Pm3m model. 
A priori, only fully randomized distribution of the pores 
is consistent with the Fm3m model (Fig. Id). Low sym- 
metry (Pm3m) becomes applicable when one finds un- 
equal occupancies for the cell body center and edge cen- 
ters (Fig. lb). These sites are abundant with the ratio 
of 1:3. If p is the occupancy of the body center, and 
1 — p/3 is the occupancy of the edge centers, the condi- 
tions p = and p = 3/4 correspond to fully ordered and 
random PB structures, respectively^. In ferrous hexa- 
cyanocobaltate, owing to different valence ratio the same 
conditions correspond to p — and p = 2/3. Intermedi- 
ate partially-ordered structures can be characterized by 
the order parameter 7 = 1 — 3p/2 . 

In Fig. 6, we have calculated the abundance of various 
coordinations of Fe 2+ for several values of 7 between 
and 1. The total abundance of the Fe 2+ species in the or- 
bital doublet ground state Ad is made up mainly of the 
abundances of FeN5H20 and irens-FeN 4 (H 2 0) 2 . Rela- 
tively small number of water-free octahedra FeN@ is as- 
sumed to be halved between two distortion isomers. 



Abundance of the orbital singlet As is com- 
posed cis-FeN 4 (H 2 0) 2 , FeN 3 (H 2 0) 3 and remaining 
FeNg (distortion isomer). Both in-plane and 3D ligand 
configurations in the 3:3 octahedron FeN 3 (H 2 0)3 con- 
tribute to A§. First, when 3 oxygens coordinating to 
iron are in one plane the x, y and z axes become all 
nonequivalent. In this case, the spin-down electron cloud 
loses the axial symmetry, and the t 2g group splits into 
three singlets. Second, when three oxygens are in 3D 
configurations, all three axes x, y and z are equivalent; 
nevertheless, the cubic symmetry is lost. Since H 2 is 
a weaker ligand than N this geometry can be viewed as 
an octahedron with three ligands expelled to longer dis- 
tance than other three. Equivalent spatial distortion is 
the elongation along the trigonal axis. Such a distortion 
was shown to stabilize the orbital singlet lowestP. 

As 7 decreases, the integral value of Ad shows in 
Fig. 6 the percentages of 63%, 52%, 44%, 39% and 
38% for 7 = 1,|,5,T,Q , respectively. Main trend of 
reduction of Ad is determined by the reduction of the 
trcms-subpopulation of FeN 4 (H 2 0) 2 . The orbital sin- 
glet ground state prevails in the disordered structure de- 
scribed by the Fm3m model. In this case (7 = 0), the ex- 
pected area of the external doublet A$ reaches the value 
of 62% similar to experimental result of 64% in Table II. 




Number of oxygens coordinating to iron 



FIG. 6: Evolution of the distribution of abundances of the 
coordination polyhedra FeNi_ n O n (n = 0,1,2,3,4,5,6) as 
the structure randomizes starting from order parameter 7 = 1 
(Pm3m model, Figl,b) to 7 = (Fm3m model, Fig.l d). Here 
p — 2(1 — 7)/3 is the occupancy of the center of the cubic cell 
in Fig.l (b). Total abundances of orbital singlet and orbital 
doublet ground states are distinguished by color. 



C. Models of short-range order in PB analogues 

The concentration of pores in archetypal Prussian 
Blue, ferric ferrocyanide Fe[Fe(CN) 6 ]i_,5-3.5H 2 0, is S = 
1/4, so that the ratio 1 — 6:6 is 3: 1. The ordered 
Pm3m structure proposed for Prussian Blue in RefP^ 
contains the high-spin cations only in coordinations FeNg 
(25%) and trans-FeN 4 (H 2 0) 2 (75%). In ferrous hcxa- 
cyanocobaltate, the ratio 1 — 6:5 is 2: 1. Therefore, 
for a vacancy-ordered phase of such a stoichiometry one 
cannot expect the same set of coordinations as for the 
Pm3m structure of the Reff^. 

In order to construct an ordered phase with the ra- 
tio 1 — 6 : (5 of 2: 1 we try to fill 3 sites with Fe, 2 
sites with Co, 1 site with Co vacancies, and sites with 
Fe vacancies. We find that among cubic structures of 
NaCl-type there is no superstructure with exact site ra- 
tio 3:2:1:0, however, a unique superstructure with the 
ratio of 30:20:12:2 exists. Using the body-centered Im3m 
model we can develop the proper arrangement of the va- 



cancies with a variable order parameter. Table III shows 
the details of the fully ordered (7 = 1) Im3m structure. 
The order parameter 7 = 1 — 8p/3 is defined via the occu- 
pancy p of the (12e) position that is the core of complex 
anion [Co(CN) 6 ] 3 ~. 

Table III. Wyckoff positions, atom type, relative co- 
ordinates and occupancies of the Na and CI nodes in 
the superstructure Im3m model on the basis NaCl-type 
framework. 



Wyckoff Atom 



2a 

6b 

8c 

12d 

12e 

24h 



Fe 



x 




y 




Occ. 




Fe Trans 1/2 1/2 1 



Co 
Co 
Co 

Fe Cis 



1/4 1/4 1/4 1 



1/4 
1/4 



1/2 




1/4 1/4 1 



Among total 14 vacancies per cubic cell of the Im3m 
structure (Table III) we place 12 vacancies of the com- 
plex anion Co(CN)g and 2 cationic (Fe 2+ ) vacancies. We 
obtained the ordered structure of Fe[Co(CN) 6 ] 2 / 3 • ffhbO 
which contains 2 very large cross-shaped pores (6+1- 
pores, or +- pores). Each cationic vacancy is surrounded 
by 6 anionic vacancies. The hydration degree || is calcu- 
lated assuming 14 water molecules per standard pore and 
an additional water molecule per Fe 2+ vacancy. Thus, 
the 6+1-pore is composed of 6 standard pores aggregated 
into 3D crosslike body. The edge of the Im3m cell is twice 
longer than the edge of the Pm3m cell and the cell volume 
is 8-fold. 

The body-centered Im3m model exhibits only two 
kinds of Fe 2+ sites with the abundance ratio of 3:1, simi- 
lar to that of the Pm3m model. In both models, the major 
site (75%) is coordinated with two waters. However, in 
the Pm3m model the major site is coordinated by waters 
in frarjs-configuration, whereas in the Im3m model the 
major site is coordinated by waters in cis-configuration. 
The minor site in Pm3m model is coordinated by six ni- 
trogens, and in Im3m model by four nitrogens and two 
waters in irons-configuration. 

It is important that we were able to construct the 
superstructure model, which strongly favors the cis- 
configuration in agreement with prevailing external dou- 
blet in Mossbauer spectra. Although the random Fm3m 
model agrees as well with the experiment, the lack of any 
short-range order between pores appears to be a severe 
limitation of this model because the vacancies bearing 
a charge interact with each other repulsively. We must 
note also that the luck of any superstructure reflections 
in diffraction patters may not be an evidence of a missing 
short-range order, but only the evidence of lacking long- 
range order. The suggested Im3m structure is a plausible 
model of the short-range order in the mixed-valence 2:3 
PB analogues. 

Although the superstructure reflections, predicted by 
either Pm3m or Im3m models, were not observed by x- 
ray diffraction, using the probabilities of local structure 
patterns observed by such a local technique as Mossbauer 
spectroscopy one can identify the proper model of short- 



range order. The Im3m model predicts the strong pre- 
dominance (3:1) of cis- over trans- population among 
the Fe 2+ environments. The model describes the situ- 
ation when concentration of elementary vacancies is so 
high that they must be in small clusters. In contrast, the 
Pm3m model describes the limit when all vacancies are 
still isolated, however, no more vacancies can be added 
without clustering. If they are actually added to pro- 
ceed from S = 1/4 to S = 1/3 (as in the bottom panel 
of Fig. 6), there will result the extra- holey Pm3m model 
that fails to describe the Mossbauer spectral data. In our 
experiments, we observe the prevailing abundance of the 
ds-configurations related to the extended crosslike pores. 
The aggregation of small pores into extended pores un- 
derlies the large (up to 750 m 2 /g) specific surface and 
unique gas sorption properties of these materials 3 ^. 



D. Charge-transferred states in the PB analogues 

The charge-transferred states which were discussed 
above to induce the distortion isomerism are seen in the 
third small (5%) doublet with IS of Fe 3 + in Fig.5 (c) 
and (e). A slight asymmetry appears also in Fig. 3, al- 
though these spectra were fitted without introducing ex- 
plicitly the Fe 3+ component. The pairs of Fe 2+ /Co 3+ 
and Fc 3+ /Co 2+ can be considered as valence tautomers, 
and we explain the asymmetry in these spectra by elec- 
tron transfer from Fe 2+ to Co 3+ . The formation and rup- 
ture of the chemical bonds implies no migration of atoms 
but only migration of electrons. The cobalt ions remain 
in low-spin state as the sixth electron of Fe 2+ transfers 
from t 2g orbital of Fe 2+ to e g orbital of Co 2+ . Other com- 
pound su scepti ble to valence tautomerism were known 
previousl} * 19 * 34 *. For example, the Prussian Blue analogue 
RbMn[Fe(CN)e]-H20, when synthesized in excess of Rb 
to fill enough framework interstitials with Rb, crystallizes 
in the structure shown in Fig. (la). In this system, the 
tautomers Fe 2+ /Mn 3+ and Fe 3+ /Mn 2+ are able to form 
different crystal phases which give way to each other as 
temperature changes. 

In the ferrous cobalticyanides, no abrupt phase change 
was found to be caused by the charge-transferred states. 
On the other hand, we observed that the content of 
Fe 3+ can be changed depending on the excess of FeCl2 
and dilution degree in the course of precipitation of 
Fe[Co(CN) 6 ] 2 /3 • 5H 2 0. Rapidly synthesized compound 
contained up to 15-20% of Fe 3+ , however, the content of 
Fe 3+ was reduced down to just 5% through reducing the 
reaction rate in dilute solutions. Most probably, the lack 
or misplacement of water molecules is needed to stabi- 
lize locally the charge-transferred states and make them 
observable in Mossbauer spectra. It appears in this case 
that the rapidly prepared samples are not saturated in C- 
water. The sites of Fe 2+ nearest to the nodes lacking the 
coordination water are subject of local tautomer trans- 
formation. In a slowly synthesized Fei.5[Co(CN) 6 ]-5H 2 
such crystallization defects are reduced. 
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Measuring the Mossbauer spectra above room temper- 
ature supports this concept. With raising temperature 
the area of the Fe 3+ subspectrum starts to increase ir- 
reversibly. Thermogravimetric analysis showed that the 
loss of sample weight occurs at heating practically start- 
ing from room temperature. The Mossbauer spectrum 
measured in situ in a furnace at 318 K shows the area of 
the Fe 3+ doublet increased up to 37%. The same sample 
measured ex-situ after prolonged (1 week) heat treatment 
at 60°C showed the area of the Fe 3+ increased up to 57%. 
The equilibrium hydration degree at 60° C is related with 
nearly half iron species transformed to Fe 3+ . Longer ex- 
posure in air at at 60°C induces no further increase of the 
spectral asymmetry and the area of spectral component 
with IS of 0.39 mm/s remains at the level of 57%. 

The isomer shift of 0.39 mm/s is attributable, in prin- 
ciple, either to high-spin Fe 3+ or to low spin Fe 2+ . 
The latter was observed below the spin-state transition 
in CsFe[Cr(CN)6p^. Theoretical approaches to lattice- 
dynamic calculations^,^ show that the change in the 
electronic state is closely related to the structure ex- 
pansion as the lattice parameter grows between the low- 
spin and high-spin forms from 1.03 nm to 1.07 nm. In 
Fe3[Co(CN)g]2 ■ nH^O, the thermal lattice expansion is 
negative^. T o veri fy that we do not deal with a reverse 
spin transitio n 1 37 ! 38 !, we have compared the magnetic sus- 
ceptibilities of the samples having 5% and 57 % of the 
subspectrum with IS of 0.39 mm/s . In the heat-treated 
sample, the susceptibility increased by 10% that is too 
small change compared to what is expected for the in- 
crease of Peff from 4.9 to 5.9 /iB- With the weight 
loss taken into account this result would rather indicate 
the formation of an intermediate-spin state of Fe 3+ which 
is known for strongly distorted octahedral coordinations 
with two alternate ligands in cis-configuratiorPn 

IV. CONCLUSION 

Mossbauer spectroscopy was applied to study the 
mixed- valence (Fe 2+ -Co 3+ ) PB analogues which dis- 
played the mixed-orbital ground states of Fe 2+ ions. The 
ground state mixing originates from both the distortion 
isomerism and the varied composition of Fe 2+ nearest- 
neighbor environment. In the A-free PB analogues, water 
molecules enter the first coordination sphere of the Fe 2+ 
ions. There appears two kinds of anisotropic ligand con- 
figurations which are the intrinsic generant of the singlet 
and doublet orbital ground states. The Fe 2+ sites are 
thus classified into two types. 

Distribution of the anionic vacancies (water-filled 
pores) over the framework is characterized by the or- 
der parameter for two long-range order models Pm3m or 



Im3m. These models are relevant to the PB analogues of 
the 4:3 and 3:2 stoichiometries, respectively. Two water 
molecules within the first coordination sphere of Fe 2+ is 
the most populated configuration of the random model 
Fm3m and ordered models Pm3m or Im3m. The short- 
range order event occurrence in PB analogues of the 3:2 
and 4:3 stoichiometries favors the cis- and trans- con- 
formations, respectively. Orbital singlet is the lowest 
in energy state for the cis-conformation. Accordingly, 
in the A-free ferrous cobalticyanide PB analogue, the 
Mossbauer spectra of Fe 2+ are dominated by the larger 
quadrupole-splitting component . 

In the A-fillcd PB analogues, where the the water lig- 
ands are withdrawn from the Fe 2+ coordination sphere, 
the least anisotropic charge distribution attainable for 
the cases of A=Rb and Cs is the orbital doublet ground 
state. For A=Na and K, the thermally-induced in- 
terconversions between the doublet and singlet ground 
states were observed attributed to distortion isomerism. 
Charge-transfered states intrinsic in the PB analogues 
with negative thermal expansion^ are suggested to be 
the source of the lattice distortions associated with the 
distortion isomerism. Asymmetry of Mossbauer spectra 
is attributed to the admixture of a small content (at the 
level of ~5%) of the Fe 3+ species associated with the 
charge-transferred states. 

Since the octahedral sites FeN 6 possess a higher sym- 
metry than the sites FeN 6 _„(H 2 0)„ we observe upon the 
insertion the A-ions into A free compounds that the value 
of Mossbauer QS drops in ferrous cobalticyanides sim- 
ilarly to the archetypal Prussian Blue Fe4[Fe(CN) 6 ] 3 • 
14H 2 0221. However, when the large-size A-ions (Cs, 
Rb) replace the small-size ones (Na,K) through ionic ex- 
change we observe the larger QS associated with lower 
symmetry local distortion isomer; the sixth d-electron of 
the Fe 2+ cation is self-localized on a single sublevel (or- 
bital singlet) of the t2 g group owing to strong lattice dis- 
tortions around the implanting sites of the low-mobility 
ions. 
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